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The behavior of quantum fluids ( 4 He and H2) within nanopores is explored in various regimes, 
using several different methods. A focus is the evolution of each fluid's behavior as pore radius R is 
increased. Results are derived with the path integral Monte Carlo method for the finite temperature 
(T) behavior of quasi-one-dimensional (ID) liquid 4 He and liquid H2, within pores of varying R. 
Results are also obtained, using a density functional method, for the T = behavior of 4 He in pores 
of variable 7?. 
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I. INTRODUCTION 

The behavior of quantum fluids (e.g. helium or hydrogen) in reduced dimensionality is one of the most interesting 
topics in low temperature physicsi Among the systems of particular interest are monolayer films (which are effectively 
two-dimensional, or 2D) and fluids within narrow pores ( ID) or droplets (OD) Properties associated with phase 

transitions (e.g. critical exponents) are particularly sensitive to the dimensionality, as exemplified by the Kosterlitz- 
Thouless theory of the 2D superfluid transition in He films. 8 Confined ideal gases are particularly amenable to 
theoretical study, because one need solve just the single particle Schrodinger equation in the given environment, rather 
than the fully interacting many-body problemJ Such an idealized model can provide qualitatively, and sometimes 
quantitatively, reliable predictions of the effects of confinement on the fluid. Furthermore, these calculations can be 
supplemented, in some cases, by virial expansions in order to assess the role of interactions, presuming their effects to 
be weaki^ii In contrast, dense quantum fluids require the use of the path integral Monte Carlo (PIMC) or diffusion 
Monte Carlo methods, which are computationally demanding approaches. Nevertheless, these techniques has been 
applied in recent years to many problems of reduced dimensionality . 1 1 1 12 

Considerable progress in understanding D < 3 behavior has also been achieved on the experimental front. Although 
the evolution of superfluidity in thin films was first observed in heat capacity experiments nearly 70 years ago, that 
subject remains incompletely understood^ The problem of quantum droplets has received particular attention in 
recent years, stimulated especially by the search for superfluid hydrogen i 6 i 14 ' 15 The discovery of carbon nanotubes 
and interest in other quasi-lD porous media have led to many investigations of ID quantum fluids^ ' 16 ? 17 

This paper addresses a set of problems involving helium and hydrogen within small radius (R) cylindrical pores, 
such as carbon nanotubes. For R ~ 0, one encounters the strictly ID limit of a quantum fluid. While there has 
been considerable exploration of the zero temperature (T) behavior of such systems, relatively few studies of finite T 
behavior of the interacting system have been carried outi 12 i 17 ' 18 i 19 i 20 i 21 i 33 In a recent study, we explored such behavior 
for ID 4 He. It was demonstrated that a ID phonon model, analogous to the Landau hydrodynamic theory of liquid 
4 He, does describe the behavior at low T of ID 4 He at moderate to high density (compared to the equilibrium density 
p = O.OSG/A)^ Section[II]of the present paper presents a finite range density functional study of the behavior of 4 He, 
at temperature T — 0, as R is varied. The results exhibit the evolution of the system from an essentially ID "axial 
phase" to a "cylindrical shell" phase, localized near the pore wallsi 16 ' 23 Section HTT1 shows how the behavior of H2 at 
finite T, but low density, varies with pore radius R. Section HVl summarizes our results. 

II. LIQUID 4 HE IN A NANOPORE AT ZERO TEMPERATURE 

As shown in previous work, finite-range density functional (FRDF) calculations permit the calculation of the 
equation of state (EOS) of helium atoms in various environments (see, e.g., Barranco et alr^ and references therein). 
In this section, we investigate the condensation of helium in a cylindrical pore of radius R and length L, with the 
same density functional employed in Refs. [2^l25ll26l with detailed description and parameters found in Ref. [27|. The 
calculation is based on the energy density £[p(r)], where p(r) is the particle density, normalized to the total number 



2 



-20 



-40 




V(0, R) 
V 

rnin 

Ho 

£ £ (Planar surface) 
V (Planar surface) 



10 
R 



] ? 



20 



FIG. 1: Cylindrical potential V^O, R) on the axis, potential V m in at the absolute minimum, and ground state energy £ gs of a 
4 He atom. For reference, the bulk chemical potential po, the minimum potential of a planar Cs sheet V^ni n (planar surface) and 
corresponding energy e gs (planar surface) are also displayed. 



of helium atoms N so that, in the pore geometry, the linear density (over 27r), n — N/(2ttL), is 

n = J dr rp(r) . (1) 

The total energy e = E/(2ttL) per unit length (over 2ir) is then 

drr£[p(r)}. (2) 

The Euler-Lagrange equation for the density is obtained by functional differentiation of this energy with respect to 
p(r), at fixed total number of particles, and takes the form 



(3) 



with fj, the chemical potential, U(p) = SE/Sp(r) the density-dependent mean field and V(r; R) the adsorption potential 
provided by the cylindrical sheet. This function was derived by considering that a He atom at radial distance r from 
the axis of the cylinder interacts pairwise with continuum sheet of atoms, of uniform areal density 9, via a Lennard- 
Jones (LJ) potential of well-depth e and hard core diameter a. The result is expressed in terms of elliptic integrals 
defined in Ref . |H 
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Here M2 n +i is an elliptic integral over angle (f>. 

Prior to determining the EOS of liquid He in pores, we consider the binding of a single He atom to an alkali 
metal cylindrical sheet. In particular, Cs substrates of diverse shapes are excellent candidates to investigate nontrivial 
aspects of the physics of wetting. While the He-Cs potential is too weakly attractive to adsorb He on a flat Cs surface 
(below saturation chemical potential /xq), there are two reasons why He does adsorb within a pore. One is that the 
surface tension cost of pore-filling is smaller than that on a planar surface (by a factor of two2£) . The other is that 
the potential can be substantially enhanced within a concave environment, either a pore or a wedge of intermediate 
openings^ 

Illustrative results are shown in Figs. [1] and [21 In Fig. [T]we plot, as functions of cylinder radius R, the potential 
V(0, R) on the cylinder axis, the potential V m i n at the absolute minimum of V(r, R), and the ground state (gs) energy 
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FIG. 2: Upper panel: potentials of a He atom inside Cs cylinders of radius for 7? values between 5 and 10 A, from left to right in 
steps of 1 A. Lower panel: corresponding probability density of one atom in its ground state. The 3D probability distributions 
have been normalized to unit maximum value. 

s gs of one atom obtained by solving the single-particle (sp) Schrodinger equation in the cylindrically symmetric 
potential. The 4 He bulk chemical potential /io = —7.15 K is displayed as a reference, as well as the minimum 
potential Vminoc = —11.3 K of a planar Cs sheet and the corresponding gs energy £ gsoo = —5.83 K. We note that 
the absolute minimum departs from the axis at a radius of about 6.8 A. On the other hand, the sp energy drops 
below the potential on the axis around R = 7.5 A, corresponding to an approximate threshold radius for the off-axis 
displacement of the maximum probability density |^(r)| 2 . This behavior is illustrated in Fig. [51 where the upper panel 
displays the potentials due to Cs sheets of varying radius R, as functions of the radial distance r, and the lower panel 
shows the probability density of the He atom. The parameters of the He-Cs interaction are e = 2.795 K, a — 5.31 A, 
and 9 = 0.38/A 2 . The large value of a reflects the extent of the hard-core repulsion for this interaction, which is 
responsible for the weak van der Waals attraction, i.e., the small value of e. 

For the smallest radii displayed, 5 and 6 A, the probability density for this "axial phase" is seen to be concentrated 
on the axis, while for the largest radii, 9 and 10 A, the particle is localized (with spread ~ 1.2 A) about 5.6 A from 
the wall; for this "cylindrical shell" phase, the latter distance asymptotes to a = 5.31 A for the hypothetical limiting 
case R = oo. In the intermediate R regime, there remains a non-negligible probability density at r = 0, seen to be 
quite sizable for R = 7 A. Note that the atom is unbound (relative to vacuum) in a 5 A tube; in fact, according to 
Fig. [TJ the gs energy is negative for radii larger than 5.2 A. Thus, a noninteracting Bose gas will be bound within 
such pores in a quasi-lD axial phase, a filled pore regime and a quasi-2D shell phase, as R increases. For even larger 
R, the usual layering structure of helium in a cavity occurs i 16 ' 26 

Next, we explore the many-body behavior by solving Eq. [3] to obtain the chemical potential /x, the total energy, 
and the density profiles p(r), as functions of the linear density n. Significant results are presented in Fig. [3J which 
shows (i and the energy per particle for two selected pore radii, namely R = 6 and 9 A, in full and dashed lines, 
respectively. For the narrower pore, a stable quasi-one-dimensional (Q1D) system with negative grand potential per 
unit length, oj = e — /x, is found for any linear density below 0.042 A -1 , above which the system becomes unstable. For 
a sufficiently wide pore, the particle distribution has moved towards the cylinder wall and the EOS displays the typical 
2D prewetting jump at n pw — 0.15 A -1 , at which point the energy per particle becomes identical to the chemical 
potential, resulting in condensation of a stable film at this or higher linear densities. This behavior can be assessed 
by examining density profiles for various values of R at a given linear density, as displayed in Fig. U for n — 0.05 A -1 . 
In this figure, the upper panel shows both the confining potential and the total one-body field U[p(r)\ + V(r;R) 
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FIG. 3: Chemical potential (black) and energy per particle (red) for liquid He in Cs pores of radii 7? = 6 and 9 A, as functions 
of linear density. The bulk chemical potential fio is shown as a dotted line. 

experienced by particles in the fluid. The resulting particle densities are plotted in the lower panel as functions of 
distance r. 

At first sight, in Fig. |4] we note the presence of the three distinct regimes of density behavior, as encountered in 
the sp pattern displayed in Fig. [5] In addition, this plot provides evidence of the effect of cohesion among the He 
atoms and its competition with the weak adhesion to the confining walls. The difference between full and dashed 
lines is the total interaction energy (per particle) of the interacting fluid, as shown in various prior studies . 24 i 25 i 26 i 27 
This total consists of a screened LJ interaction plus two density-dependent correlation terms, one attractive, the other 
repulsive, representing the effects of the many-body environment on the one-body field. For the two smallest radii, 
the overall outcome, as seen in the upper panel of Fig. [4] is a visible reduction in potential energy and a distortion 
of the potential well that causes an earlier departure of the minimum off-axis, as compared with the sp picture in 
Fig. [21 The situation is reversed for the largest radii. For R = 9 A, the potential energy gain of about 2 K near the 
cylinder axis is sufficient to sustain a non-negligible amount of fluid in that region. As a consequence, pores of radii 
between somewhat less than 7 A and more than 9 A result in a particle density spread out from axis to wall, and 
the appearance of the purely axial regime shifts towards larger radii. For radii near 12 A and above, instead, the 
correlation effects in the helium fluid are negligible and the system behaves essentially as a confined gas. 

III. H 2 IN A NANOPORE AT FINITE T 

We consider in this section the adsorption of H2 (which, after 4 He, is the textbook quantum fluid) molecules inside 
a cylindrical pore within a bulk solid. As in the preceding section, a low density gas of molecules within a large pore 
will typically adsorb first onto the inner surface of the pore, creating a 2D cylindrical shell phase at low densities; 
at higher densities, central regions of the pore will fill with a 3D bulk adsorbate^ Very narrow pores, however, will 
constrain the adsorbate to lie near the pore axis, forming a quasi-lD ("axial") phase. This variation of behavior 
provides the opportunity to control the effective dimensionality of the adsorbate by adjusting the size of the pore into 
which it adsorbs. For concreteness, we chose to study pores in MgO glass, a much more strongly adsorbing substrate 
than the Cs discussed in the previous section; the well-depth (48 meV) of an H2 molecule on planar MgO is about 15 
times higher than that on Csi 29 i 30 MgO nanopores can be fabricated as small as 3 nm in diameter^ The H 2 -MgO 
potential was taken to be derived from Lennard- Jones 12-6 pair interactions, integrated over the bulk MgO solid; 
these gas-surface interaction parameters are given by a = 2.014 A, e = 45.91 K. The H2-H2 interaction potential 
was taken to be of the Silvera-Goldman formi^l The molecule is treated as a spherically symmetric point particle, 
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neglecting internal degrees of freedom. 



A. Structure: path integral Monte Carlo simulations 

This quantum statistical system can be simulated by the path integral Monte Carlo (PIMC) method i 11 ^ 2 The 
simulations were performed for various pore radii R at fixed temperature T = 0.5 K (or 1 K) and particle number 
N = 28; the mean 3D density p = N/V changed as the pore volume was varied. The value N — 28 was chosen to 
approximate the equilibrium linear density of a ID system of fixed length (of the periodic unit cell) L = 125 A . 33 > 34 
The value of L was chosen large enough to ensure that finite size effects were minimal, as confirmed by simulation 
tests.— 

The distribution of H2 molecules calculated by PIMC simulation is given in Fig. [5] for several pore sizes. The 
axial-to-shell transition is more apparent in Figure [6J which shows the H2 radial density distribution, along with the 
corresponding H2-pore potential energy curves. Classically, the molecules will concentrate on axis when the potential 
minimum is at r = 0, and will move off axis when the minimum is found at some r > 0, if intermolecular interactions 
are neglected (the low density limit). The PIMC simulations at low densities indicate that the adsorbate moves off 
axis between pore radii of R = 2.5 and 2.75 A; by comparison, the classical potential minimum moves off axis at 
slightly smaller pore size, R — 2.35 A. In other words, the molecules tend to concentrate near the axis more when 
treated as quantum particles than when treated classically. The comparison between quantum and classical behavior 
in this system is discussed more fully in the following section. 

A closer examination of samples drawn from Monte Carlo simulation, Fig. [JJ reveals the presence of bound H2-H2 
pairs, or "dimers" , as well as possible bound triplets, or "trimers" . Such clustering is expected below a temperature 
set by the H2-H2 binding energy. This visual evidence is corroborated quantitatively by the pair correlation function 
derived from PIMC simulation, Fig. [8l which shows an enhanced probability for H2 molecules to be separated by 
r = 4.8 A. By comparison, the classical equilibrium separation of a pair of H2 molecules is r — 3.4 A; the separation 
of quantum H2 is greater due to zero-point fluctuations, which reduce the binding energy to about 3 K, roughly 
one-tenth of the well-depth of the pair potential i 31 ^ 36 

In a crude analysis of Fig. [71 we estimate that at T = 1 K approximately equal numbers of H2 molecules participate 
in monomers (unpaired H2) as dimers, i.e. a dimer/monomer ratio of 0.5. This quantity can be roughly estimated by 
assuming that the monomer<^>dimer process is in detailed balance, implying that their respective chemical potentials, 
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FIG. 5: H2 density distribution [at T — 0.5 K, p = (2.55 x 1CP S A)/R 2 ], viewed in transverse section, for pore radii 7? = 2.5- 
4.5 A, as indicated. Depicted are several hundred superimposed Monte Carlo samples of particle positions (dots) for a system 
of N = 28 particles, along with the inner surface of the pore (circles). (The circles are distorted into ellipses due to the aspect 
ratio of the figure.) The H2 evolves from a ID (axial) to a 2D (cylindrical shell) phase as R increases. 



Hi and fi2, satisfy {12 = 2 fix- Assuming the H2 to be an ideal gas of N molecules, the chemical potentials are related 
to the individual populations (JV») by 



fix=k B T\n[— — j, 



M2 = k B Th\ 



eV 
eV 



(5) 



where A = y 2irh 2 / (rnksT) is the thermal wavelength and e is the binding energy per dimer. These relations, under 
the constraint N — N% + 2N2, give a dimer/monomer ratio 



7V 2 1 + 4 X - ^/i + 8x 



2 D/2f±\ e e/(k B T) ^ (g) 



Nx 2(vT+8jb-1) Vn Q 

Here n — N/V is the mean 3D number density and jiq = 1/X D is the quantum concentration. For N — 28 H2 
molecules at T = 1 K in a 3D volume V = ttR 2 L, where R = 7 A and L = 125 A, assuming their binding energy is 
e = 3 K, Eq. [6] predicts a dimer-monomer ratio of 8.5. This estimate exceeds the ratio observed in PIMC simulation by 
an order of magnitude. However, the preceding back-of-the-envelope calculation is expected to have limited accuracy 
since it ignores the quasi-2D geometry of adsorption, the presence of internal rotational- vibrational excitational modes 
of the dimer and the possibility of trimer, or larger, clusters of bound molecules. 



B. Structure: semiclassical approximation 



As mentioned above, the axial and shell phases can be distinguished in the classical, noninteracting case by the 
location of the potential minimum (r = or r > 0) for a given pore radius R. Quantum mechanically, one must 




FIG. 6: H 2 probability density [at T = 0.5 K, p = (2.55 x 1(T 3 k)/R 2 } and MgO pore potential, for pore radii R = 2.5, 2.75, 
3.0 and 3.25 A, as functions of dimensionless radius r/R. (Radial densities near r — are exaggerated due to finite size effects 
after normalizing the radial distribution by l/(27rr) to obtain the probability density.) 
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FIG. 7: H2 spatial density distribution (at T = 1 K, average density p — 1.46 x 10 _3 /A 3 ), for N = 28 particles in a pore 
of radius 7 A. Depicted are two individual Monte Carlo samples of particle positions. (The scale of the vertical, longitudinal 
axis is compressed relative to that of the transverse axes.) Note the presence of dimers (H2-H2 pairs) in both plots, as well as 
possible trimers. 



resort to either solving the Schrodinger equation and finding the peak probability from the radial wavefunction (for 
the noninteracting case) or by carrying out a full path integral Monte Carlo simulation (in the general interacting 
case). 

Besides their use in Monte Carlo simulations, the path integral method also yields a semiclassical approximation to 
the density that is almost as simple as inspection of the classical potential. A variational approximation, originally due 
to Feynman, reduces the quantum statistical path integral to the ordinary Boltzmann integral of classical statistical 
mechanics, with the classical potential V replaced by a semiclassical effective potential V c g containing first-order 
quantum corrections to VM- This potential is expressed as a weighted average of the real potential over a distance 
comparable to the de Broglie thermal wavelength. Thus, in this approximation, the axial phase of quantum H2 occurs 
when the minimum of the effective potential lies at r = 0. Feynman's variational approximation yields a formula for 
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FIG. 8: H2 pair correlation function (at T = 1 K, p = 1.46 x 10 3 /A 3 ), for N = 28 particles in a pore of radius 7 A. Inset: 
expanded view of the pair correlation function, near the r = 4.8 A peak, that indicates the presence of dimerization. 



An 



V e ff as a functional of the classical potential in the case of a noninteracting quantum gas, 

1 



V eS (r) = (2^A 2 )- 3/2 / d 3 r' V(r') exp [--(|r - r'|/A) 2 J , A = h/y/l2mk B T. 
Here m is the mass of a particle, in this case an H2 molecule. This has a high-temperature expansion 

V eS = V + ^A 2 V 2 y + 0(A 4 ) . 



(7) 



(8) 



The effective potential in Eq. [7] is the convolution of the classical potential with a Gaussian of width A, and may 
be interpreted as a "quantum smeared" version of the classical potential. In the high temperature (small A) limit, 
the smearing vanishes, and the effective potential reproduces the classical result. In the case of cylindrical symmetry, 
these equations become 



V eB {r) = A" 



! ^°°dr'y(rV/o(^)exp 



-i(r 2 +r' 2 )/A 2 



(9) 



and 



V c ff = V + -A 2 (V" + V'/r) + 0(A 4 ) 



(10) 



Effective potentials for various R and T values are displayed in Figure[3J As expected, V e g approaches V as T increases. 
The main qualitative results concern the location of the effective potential minimum. This energy is greater than 
that of the classical potential minimum, due to the zero point energy. It is also located closer to the axis, because the 
repulsive hard wall of the potential at r = R damps the wavefunctions of quantum particles in the repulsive region. 

Qualitatively, then, these semiclassical considerations imply that a quantum gas has a greater tendency to be 
confined near the axis than does a classical gas, agreeing with the preceding results obtained by PIMC simulation. 
Quantitatively, the semiclassical approximation breaks down at temperatures above the T = 0.5 K studied by PIMC, 
but one may use it for qualitative purposes. For example, the location of the effective potential minimum at higher 
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FIG. 9: Hb-MgO pore radial effective potentials at T = 100 and 200 K, compared to the corresponding classical potentials, at 
pore radii R = 2.2, 2.4, and 2.6 A. Note that the quantum behavior described by the effective potential approaches classical 
behavior in the high temperature limit. 



temperatures (T > 50 K) is given in Figure [TUl Note that the classical potential minimum moves off axis near 
R = 2.30-2.35 A, while the effective potential minimum moves off axis near 2.35 A at 100 K and near 2.45 A at 50 K. 
By extrapolation of these semiclassical predictions, it is plausible that the low T axial-shell transition occurs between 
2.50 and 2.75 A, as found in the preceding section by PIMC simulation at T — 0.5 K. 



C. Heat capacity 

In accordance with the structural results of the preceding sections, the heat capacity C = (dE/dT)M,v of the H2 
is expected to exhibit ID behavior for small R, and 2D behavior for larger R pores. In the ideal gas approximation, 
for small R, the single longitudinal degree of freedom should give a dimensionless specific heat, or heat capacity per 
particle (in ks = 1 units) of 1/2, in accordance with the equipartition theorem. An additional azimuthal degree of 
freedom is available in the cylindrical shell phase of larger R, increasing C /(Nks) to 1, given enough thermal energy 
to overcome the azimuthal excitation energy gap. At very high T, the large radial excitation energy gap will be 
exceeded, allowing an additional degree of freedom, raising the specific heat to 3/2. 

In the noninteracting case, the heat capacity can be calculated by direct differentiation of the average energy, 
upon obtaining the energy spectrum from numerical solution of the Schrodinger equation. In the grand canonical 
ensemble, the average particle number and energy of a Bose gas are given by summing the Bose-Einstein distribution 
of states over all possible energies (weighted by the density of states, or relative prevalence of each energy); this 
includes an integration over the continuum of longitudinal energies as well as a summation over the discrete spectrum 
of radial/azimuthal transverse energies e nyV . These quantities are then given by 
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FIG. 10: Minimum of the Hh-MgO pore radial effective potential, r m i n , as a function of pore radius R, for various temperatures 
(increasing from lower right to upper left), compared to the minimum of the classical potential (upper leftmost curve). The 
axial/shell transition occurs at the radius R for which r m i n becomes nonzero. This transition point for the quantum gas is the 
same as the classical transition of R = 2.3 A when T — 200 K, and increases to 2.45 A as the temperature of the quantum gas 
decreases to T = 50 K. (The resolution of the plot and therefore the precision of the estimated transition point is calculated to 
within 0.05 A.) 

where z is the fugacity [exp(f3E)), fi is the chemical potential, (3 — l/(ksT) is the inverse temperature, g(E) = 
{L/ \)[f3/ (itE)] 1 / 2 is the density of states for a ID (longitudinal) gas, and the degeneracy g v is 1 for v — (no angular 
momentum) and 2 for v > [corresponding to the ±(counter-)clockwise azimuthal states]. 

Fixing the particle number N, the first equation is used to solve [i in terms of particle number, fJ>(N), which is 
substituted into the second equation to determine the average energy E(T,N). The specific heat c = C/N, or heat 
capacity per particle, is obtained by a finite difference approximation to the derivative [d(E / 'N) / '8t]n,v ■ 

The resulting dimensionless specific heat as a function of R, at T = 1 and 5 K, is depicted in Figure QT] (Also 
shown are predictions for a quasi-lD approximation, in which only the azimuthal ground and first excited states are 
considered, not discussed further here.) The evolution of the specific heat from 1/2 for small R to 1 for large R is 
readily apparent. The small R specific heat is less than the classical prediction of 1/2 because the temperatures are 
low enough that the quantum system has not reached the classical limit; it is closer to 1/2 for the more classical 
T = 5 K case than for the lower temperature, T = 1 K, case. 

The higher T case approaches the 2D limit more quickly, since more thermal energy is available to excite the higher 
azimuthal modes that contribute the additional degree of freedom necessary for 2D motion. Structurally, in the PIMC 
simulation at T = 0.5 K, the H2 changed from an axial to a shell-like configuration for R < 3 A. The heat capacity, 
on the other hand, indicates a transition from ID to 2D behavior at larger 7?~6AatT=lK. Thus, the H2 may 
behave thermodynamically as a ID gas even when it is structurally arranged as a 2D film. At pore sizes of R — 4 A, 
the pore is large enough for a shell configuration to be energetically favorable, yet the azimuthal energy gap, which 
shrinks with increasing R, is not yet large enough to permit appreciable rotational motion leading to an enhanced 
heat capacity. (The energy gaps of the first few azimuthal excited states are also depicted in Fig. [TT1) 
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FIG. 11: Above: Specific heat of noninteracting H2 gas in a MgO pore as a function of pore radius R, at fixed mean linear 
density N/L = 0.224 A -1 [mean 3D density p — N/L/(tyR 2 )], for T = 1 and 5 K (right and left solid curves). Contrasted are the 
predictions of the quasi-lD model with an energy gap e = eo.i — £0,0 determined by the ground state and first azimuthal state 
(dashed curves). (The quantum numbers n, v of £„,„ denote radial and azimuthal excitations, respectively.) The horizontal 
dashed lines are the classical ID and 2D limits of 5 and 1, respectively. Below: The energy gap between the ground state 
(eo,o) an d the first four azimuthal excitations of the ground radial state (eo,{i...4}> increasing from left to right). The horizontal 
dashed lines correspond to energy gaps of 1 and 5 K. 

IV. SUMMARY 

Results have been obtained for the T = properties of liquid 4 He and finite T properties of liquid H2. In both 
cases, the behavior is seen to exhibit dimensional crossover as a function of R, as shown in the density profiles 
and thermodynamic properties (energies, chemical potentials and specific heat). A wide variety of other properties, 
including detailed dependences on N and T, are amenable to further theoretical investigation. While we hope to 
undertake some of these investigations in the future, they require a substantial increase in computational resources 
and time. Nevertheless, these are justified by recent experimental studies of quantum fluids in and near carbon 
nanotubes and other porous mediai 38 ' 39 ! 40 ' 41 
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